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I. OUTLINE 

Quantum cellular automata (QCA) are a generalization of (classical) cellular automata (CA) and in particular of 
reversible CA. The latter are reviewed shortly. An overview is given over early attempts by various authors to define 
one-dimensional QCA. These turned out to have serious shortcomings which are discussed as well. Various proposals 
subsequently put forward by a number of authors for a general definition of one- and higher-dimensional QCA are 
reviewed and their properties such as universality and reversibility are discussed. 
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II. GLOSSARY 

Space-homogeneous : The transition function / update table is the same for each cell. 
Time- homogeneous : The transition function / update table is time-independent. 
Configuration : The state of all cells at a given point in time. 

Neighborhood : All cells with respect to a given cell that can affect this cell's state at the next time step. A 
neighborhood always contains a finite number of cells. 

Update table : Takes the current state of a cell and its neighborhood as an argument and returns the cell's state 
at the next time step. 

Schrodinger picture : Time evolution is represented by a quantum state evolving in time according to a time- 
independent unitary operator acting on it. 

Heisenberg picture : Time evolution is represented by observables (elements of an operator algebra) evolving in 
time according to a unitary operator acting on them. 

BQP complexity class : Bounded error, quantum probabilistic, the class of decision problems solvable by a quantum 
computer in polynomial time with an error probability of at most 1/3. 

QMA complexity class : Quantum Merlin- Arthur, the class of decision problems such that a "yes" answer can be 
verified by a 1-message quantum interactive proof (verifiable in BQP). 

Quantum Turing machine : A quantum version of a Turing machine - an abstract computational model able to 
compute any computable sequence. 

Swap operation : The one-qubit unitary gate U = 
Hadamard gate : The one-qubit unitary gate U — -j= 
Phase gate : The one-qubit unitary gate U = 

Pauli operator : The three Pauli operators are a x = I " * ) , a y = I " n ' ) , a z = 







III. DEFINITION OF THE SUBJECT AND ITS IMPORTANCE 

Quantum cellular automata (QCA) are a quantization of classical cellular automata (CA), d-dimensional arrays of 
cells with a finite-dimensional state space and a local, spatially-homogeneous, discrete-time update rule. For QCA 
each cell is a finite-dimensional quantum system and the update rule is unitary. CA as well as some versions of 
QCA have been shown to be computationally universal. Apart from a theoretical interest in a quantized version of 
CA, QCA are a natural framework for what is most likely going to be the first application of quantum computers 
- the simulation of quantum physical systems. In particular, QCA are capable of simulating quantum dynamical 
systems whose dynamics are uncomputable by classical means. QCA are now considered one of the standard models 
of quantum computation next to quantum circuits and various types of measurement-based quantum computational 
models 1 . Unlike their classical counterpart, an axiomatic, all-encompassing definition of (higher-dimensional) QCA 
is still missing. 



IV. INTRODUCTION 

Automata theory is the study of abstract computing devices and the class of functions they can perform on their 
inputs. The original concept of cellular automata is most strongly associated with John von Neumann (*1903, +1957), 
a Hungarian mathematician who made major contributions to a vast range of fields including quantum mechanics, 
computer science, functional analysis and many others. According to Burks, an assistant of von Neumann, ((l|) von 
Neumann had posed the fundamental questions: "What kind of logical organization is sufficient for an automaton 
to reproduce itself?" . It was Stanislav Ulam who suggested to use the framework of cellular automata to answer 
this question. In 1966 von Neumann presented a detailed analysis of the above question in his book Theory of 
Self -Reproducing Automata |l|). 



For details on these and other aspects of quantum computation see the article by Kendon in this issue. 
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Thus, von Neumann initiated the field of cellular automata. He also made central contributions to the mathematical 
foundations of quantum mechanics and, in a sense von Neumann's quantum logic ideas were an early attempt at 
defining a computational model of physics. But he did not pursue this, and did not go in the directions that have led 
to modern ideas of quantum computing in general or quantum cellular automata in particular. 

The idea of quantum computation is generally attributed to Feynman who, in his now famous lecture in 1981, pro- 
posed a computational scheme based on quantum mechanical laws (0). A contemporary paper by Benioff contains the 
first proposal of a quantum Turing machine (0). The general idea was to devise a computational device based on and 
exploiting quantum phenomena that would outperform any classical computational device. These first proposals were 
sequentially operating quantum mechanical machines imitating the logical operations of classical digital computation. 
The idea of parallelizing the operations was found in classical cellular automata. However, how to translate cellular 
automata into a quantum mechanical framework turned out not to be trivial. And to a certain extent how to do this 
in general remains an open question until today. 

The study of quantum cellular automata (QCA) started with the work of Grossing and Zeilinger who coined the 
term QCA and provided a first definition (0). Watrous developed a different model of QCA (0). His work lead 
to further studies by several groups (0; 0; S). Independently of this, Margolus developed a parallelizable quantum 
computational architecture building on Feynman's original ideas (0) . For various reasons to be discussed below, none 
of these early proposals turned out to be physical. The stud y o f QCA gained new momentum with the work by 
Richter, Schumacher, and Werner (flol ; [Til ) and others (fl2l ; Il3l Il4) who avoided unphysical behavior allowed by the 
early proposals (jlll ; Il4 ) . It is important to notice that in spite of the over two-decade long history of QCA there is no 
single agreed-upon definition of QCA, in particular of higher-dimensional QCA. Nevertheless, many useful properties 
have been shown for the various models. Most importantly, quite a few models are shown to be com put ationally 
universal, i.e. they can simulate any quantum Turing machine and any quantum circuit efficiently ((Bl; \dt. Il2t HBk Il6h. 
Very recently, their ability to generate and transport entanglement has been illustrated (fl7l). 

A comment is in order on a class of models which is often labeled as QCA but in fact are classical cellular automata 
implemented in quantum mechanical structures. They do not exploit quantum effects for the actual computation. To 
make this distinction clear they are now called quantum-dot QCA. These types of QCA will not be discussed here. 



V. CELLULAR AUTOMATA 

Definition (Cellular Automata). A cellular automaton (CA) is a J^-tuple (L, £,A/\ /) consisting of (1) a d- 
dimensional lattice of cells L indexed i £ Z d , (2) a finite set of states £, (3) a finite neighborhood scheme M C 
and (4) a local transition function f : S - ^ — + S. 

A CA is discrete in time and space. It is space and time homogeneous if at each time step the same transition 
function, or update rule, is applied simultaneously to all cells. The update rule is local if for a given lattice L and 
lattice site x, f(x) is localized in x + J\f = {x + n\x G L, n G A/"}, where J\f is the neighborhood scheme of the CA. In 
addition to the locality constraint the local transition function / must generate a unique global transition function 
mapping a lattice configuration Ct G S L at time t to a new configuration Ct+i at time f + 1: F : S L — > E L . Most CA 
are defined on infinite lattices or, alternatively, on finite lattices with periodic boundary conditions. For finite CA 
only a finite number of cells is not in a quiescent state, i.e. a state that is not effected by the update. 

The most studied CA are the so-called elementary CA - 1-dimensional lattices with a set of two states and a 
neighborhood scheme of radius 1 (nearest-neighbor interaction), i.e. the state of a cell at point x at time t + 1 only 
depends on the state of cells x — 1, x, and x + 1 at time t. There are 256 such elementary CA, easily enumerated 
using a scheme invented by Wolfram Ijl8l ). As an example and for later reference, the update table of rule 110 is given 
in Table |TJ CA with update rule '110' have been shown to be computationally universal, i.e. they can simulate any 
Turing machine in polynomial time (fl9h . 

110 _ 111 110 101 100 011 010 001 000 
""01101110 

TABLE I Update table for CA rule '110' (the second row is the decimal number '110' in binary notation). 

A possible approach to constructing a QCA would be to simply "quantize" a CA by rendering the update rule 
unitary. There are two problems with this approach. One is, that applying the same unitary to each cell does not 
yield a well-defined global transition function nor necessarily a unitary one. The second problem is the synchronous 
update of all cells. "In practice", the synchronous update of, say, an elementary CA can be achieved by storing the 
current configuration in a temporary register, then update all cells with odd index in the original CA, update all cells 
with even index in the register and finally splice the updated cells together to obtain the original CA at the next 



FIG. 1 Even (solid lines) and odd (dashed lines) of a Margolus partitioning scheme in d — 2 dimensions using blocks of size 
2x2. For each partition one block is shown shaded. Update rules are applied alternatingly to the solid and dashed partition. 

time step. Quantum states, however, cannot be copied in general due to the so-called no-cloning theorem (poh . Thus, 
parallel update of a QCA in this way is not possible. Sequential update on the other hand leads to either an infinite 
number of time steps for each update or inconsistencies at the boundaries. One solution is a partitioning scheme as 
it is used in the construction of reversible CA. 



A. Reversible Cellular Automata 

Definition (Reversible CA) . A CA is said to be reversible if for every current configuration there is exactly one 
previous configuration. 

The global transition function F of a reversible CA is bijective. In general, CA are not reversible. Only 16 out of 
the 256 elementary CA rules are reversible. However, one can construct a reversible CA using a partitioning scheme 
developed by Toffoli and Margolus for 2-dimensional CA (|2lh . 

Consider a 2-dimensional CA with nearest neighborhood scheme Af = {x € Z 2 |V|a;i| < 1}. In the partitioning 
scheme introduced by Toffoli and Margolus each block of 2 x 2 cells forms a unit cube □ such that the even translates 
□ + 2x with x £ Z 2 and the odd translates □ + 1 + 2x, respectively, form a partition of the lattice, see Fig. [T] 
The update rule of a partitioned CA takes as input an entire block of cells and outputs the updated state of the 
entire block. The rule is then applied alternatingly to the even and to the odd translates. The Margolus partitioning 
scheme is easily extended to d-dimensional lattices. A generalized Margolus scheme was introduced by Schumacher 
and Werner (fill). It allows for different cell sizes in the intermediate step. 

A partitioned CA is then a CA with a partitioning scheme such that the set of cells are partitioned in some periodic 
way: Every cell belongs to exactly one block, and any two blocks are connected by a lattice translation. Such a CA is 
neither time homogeneous nor space homogeneous anymore, but periodic in time and space. As long as the rule for 
evolving each block is reversible, the entire automaton will be reversible. 

VI. EARLY PROPOSALS 

Grossing and Zeilinger were the first to coin the term and formalize a QCA Qj. In the Schrodinger picture of 
quantum mechanics the state of a system at some time t is described by a state vector \ip t ) in Hilbert space H. The 
state vector evolves unitarily, 

\^ t+1 ) = U\^ t ) . (1) 

U is a unitary operator, i.e. UW = 1, with the complex conjugate IP and the identity matrix 1. If {|<fo)} is a 
computational basis of the Hilbert space TL any state € tL can be written as a superposition Y^ua c i \4>i)i with 
coefficients Cj € C and X^i c i c l = 1- The QCA constructed by Grossing and Zeilinger is an infinite 1-dimensional 
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lattice where at time t lattice site i is assigned the complex amplitude Ci of state \ipt}- The update rule is given by 
unitary operator U. 

Definition (Grossing-Zeilinger QCA). A Grossing-Zeilinger QCA is a 3-tuple (L,H,U) which consists of (1) an 
infinite 1- dimensional lattice L C Z representing basis states of (2) a Hilbert space H with basis set {\<f>i)}, and (3) a 
band- diagonal unitary operator U . 

Band-diagonality of U corresponds to a locality condition. It turns out that there is no Grossing-Zeilinger QCA with 
nearest-neighbor interaction and nontrivial dynamics. In fact, later on, Meyer showed more generally that "in one 
dimension there exists no nontrivial homogeneous, local, scalar QCA. More explicitly, every band r-diagonal unitary 
matrix U which commutes with the one-step translation matrix T is also a translation matrix T k for some k G Z, 
times a phase" (f22l ). 

Grossing and Zeilinger also introduce QCA where the unitarity constraint is relaxed to only approximate unitarity. 
After each update the configuration can be normalized which effectively causes non-local interactions. 

The properties of Grossing-Zeilinger QCA are studied by Grossing and co-workers in some more detail in following 
years, see (|23l ) and references therein. This pioneering definition of QCA, however, has not been studied much 
further, mostly because the "non-local" behavior renders the Grossing-Zeilinger definition non-physical. In addition, 
it has little in common with the concepts developed in quantum computation later on. The Grossing-Zeilinger 
definition really concerns what one would call today a quantum random walk (for further details see the review by 
Kempe (H3)). 

The first model of QCA researched in depth was that introduced b y W atrous (0), whose ideas where further explored 
by van Dam (@), Diirr, LeThanh, and Santha (0; H), and Arrighi l|25tl . A Watrous-QCA is defined over an infinite 
1-dimensional lattice, a finite set of states including a quiescent state. The transition function maps a neighborhood 
of cells to a single quantum state instantaneously and simultaneously. 

Definition (Watrous-QCA). A Watrous-QCA is a 4-tuple (h,T,,Af,f) which consists of (1) a 1-dimensional lattice 
L C Z, (2) a finite set of cell states E including a quiescent state e, (3) a finite neighborhood scheme N ' , and (4) a 
local transition function f : E^ — > TL^, ■ 

Here, 7is denotes the Hilbert space spanned by the cell states E. This model can be viewed as a direct quantization 
of a CA where the set of possible configurations of the C A is extended to include all linear superpositions of the classical 
cell configurations, and the local transition function now maps the cell configurations of a given neighborhood to a 
quantum state. One cell is labeled "accept" cell. The quiescent state is used to allow only a finite number of states to 
be active and renders the lattice effectively finite. This is crucial to avoid an infinite product of unitaries and, thus, 
to obtain a well-defined QCA. 

The Watrous QCA, however, allows for non-physical dynamics. It is possible to define transition functions that do 
not represent unitary evolution of the configuration, either by producing superpositions of configurations which do 
not preserve the norm, or by inducing a global transition function which is not unitary. This leads to non-physical 
pro perties such as super-luminal signaling (fill ) . The set of Watrous QCA is not closed under composition and inverse 

Watrous defines a restricted class of QCA by introducing a partitioning scheme. 

Definition (Partitioned Watrous QCA). A partitioned Watrous QCA is a Watrous QCA with E = E| x E c x S r for 

finite sets E;, E c , and E r , and matrix A of size E x E. For any state s = (s;, s c , s r ) 6 E define transition function f 
as 

f(si,s 2 ,s 3 ,s) = A( Si3iSm2jSri)iS , (2) 

with matrix element A SiiS .. 

In a partitioned Watrous QCA each cell is divided into three sub-cells - left, center, and right. The neighborhood 
scheme is then a nearest-neighbor interaction confined to each cell. The transition function consists of a unitary acting 
on each partitioned cell and swap operations among sub-cells of different cells. Fig. [2] illustrates the swap operation 
between neighboring cells. 

For the class of partitioned Watrous QCA Watrous provides the first proof of computational universality of a 
QCA by showing that any quantum Turing machine can be efficiently simulated by a partitioned Watrous-QCA with 
constant slowdown and that any partitioned Watrous-QCA can be simulated by a quantum Turing machine with 
linear slowdown. 

Theorem 1 ((d)). Given any quantum Turing machine Mtm, there exists a partitioned Watrous QCA Mqa which 
simulates Mtm with constant slowdown. 
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FIG. 2 Each cell is divided into three sub-cells labeled 1, c, and r for left, center, and a right, respectively. The update rule 
consists of swapping left and right sub-cells of neighboring cells and then updating each cell internally using a unitary operation 
acting on the left, center, and right part of each cell. 

Theorem 2 ((|5|)). Given any partitioned Watrous QCA Mqa, there exists a quantum Turing machine Mtm which 
simulates Mqa with linear slowdown. 

Watrous' model was further developed by van Dam (@), who defines a QCA as an assignment of a product vector 
to every basis state in the computational basis. Here the quiescent state is eliminated and thus the QCA is made 
explicitly finite. Van Dam showed that the finite version is also computationally universal. Efficient algorithms to 
decide whether a given 1-dimensional QCA is unitary was presented by Durr, LeTanh, and Santha Due to 

substantial shortcomings such as non-physical behavior, these early proposals were replaced by a second wave of 
proposals to be discussed below. 

VII. DEFINITION 

There is not a generally accepted QCA model that has all the attributes of the CA model: unique definition, simple 
to describe, and computationally powerful. In particular, there is no axiomatic definition, contrary to its classical 
counterpart, that yields an immediate way of constructing/enumerating all of the instances of this model. Rather, 
each set of authors defines QCA in their own particular fashion. In general, a QCA consists of a d-dimensional lattice 
of identical finite-dimensional quantum systems, a finite set of states, a finite neighborhood scheme, and a set of local 
unitary transition rules. The states s£E are basis states spanning a finite-dimensional Hilbert space. At each point 
in time a cell represents a finite-dimensional quantum system in a superposition of basis states. The unitary operators 
represent the discrete-time evolution of strictly finite propagation speed. 

VIII. MODELS OF QCA 
A. Reversible QCA 

Schumacher and Werner use the Hcisenberg picture rather than the Schrodinger picture in their model (fill). Thus, 
instead of associating a d-level quantum system with each cell they associate an observable algebra with each cell. 
Taking a quasi-local algebra as the tensor product of observable algebras over a finite subset of cells, a QCA is then 
a homomorphism of the quasi-local algebra, which commutes with lattice translations and satisfies locality on the 
neighborhood. 

The observable-based approach was first used in Ref. (| 10h with focus on the irreversible case. However, this definition 
left questions open such as whether the composition of two QCA will again form a QCA. The following definition 
does avoid this uncertainty. 

Consider an infinite d-dimensional lattice L C Z d of cells x € Z d , where each cell is associated with the observable 
algebra A x and each of these algebras is an isomorphic copy of the algebra of complex d x d-matrices. When A C It 
is a finite subset of cells, denote by .4(A) the algebra of observables belonging to all cells in A, i.e. the tensor product 
(E>x€aA x . The completion of this algebra is called a quasi-local algebra and will be denoted by ACL d ). 

Definition (Reversible QCA). A quantum cellular automaton with neighborhood scheme M C II 1 is a homomorphism 
T : A(1* d ) — > A{1, d ) of the quasi-local algebra, which commutes with lattice translations, and satisfies the locality 
condition T(A(A)) C T(A(A + AT)) for every finite set A C Z d . The local transition rule of a cellular automaton is 
the homomorphism Tq : Aq — > A(J\f). 

They present and prove the following theorem on one-dimensional QCA. 
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FIG. 3 Generalized Margolus partitioning scheme in 1 dimension using two unitary operations U and V. 



Theorem 3 (Structure Theorem (fllh ). Let T be the global transition homomorphism of a one-dimensional nearest- 
neighbor QCA on the lattice Z d with single-cell algebra Aq = M.d- Then T can be represented in the generalized 
Margolus partitioning scheme, i.e. T restricts to an isomorphism 

T:4D)^0B 8 , (3) 

where for each quadrant vector q £ Q, the subalgebra B q C A(0 + q) is a full matrix algebra, B q A4 n ^ . These algrebras 
and the matrix dimensions n(q) are uniquely determined by T. 

Theorem [3] does not hold in higher dimensions (f26h . 

A central result obtained in this framework is that almost any |26l ) 1-dimensional QCA can be represented using a 
set of local unitary operators and a generalized Margolus partitioning (jllh. as illustrated in Fig. [3] Furthermore, if the 
local implementation allows local ancillas, then any QCA, in any lattice dimension can be built from local unitaries 
((Til; l2ot). In addition, they prove the following Corollary. 

Corollary 1 ((11)). The inverse of a nearest-neighbor QCA exists, and is a nearest-neighbor QCA. 

The latter result is not true for CA. A similar result for finite configurations is obtained in (fl4h . Here evidence 
is presented that the result does not hold for two dimensional QCA. The work by Schumacher and Werner can be 
considered the first general definition for 1-dimensional QCA. A similar result for many-dimensional QCA does not 
exist. 



B. Local unitary QCA 

Peres-Delgado and Cheung propose a local unitary QCA (12). 

Definition (Local-unitary QCA). A local-unitary QCA is a 5-tuple {(L, S, N, C/q, Vo)} consisting of (1) a d- 
dimensional lattice of cells indexed by integer tuples L C Z, d , (2) a finite set of orthogonal basis states E, (3) a 
finite neighborhood scheme AT C Z d , (4) a local read function Uq : {TL^,)®^ — * (We)®"^, and (5) a local update 
function Vo : 7is — * TLy,- The read operation carries the further restriction that any two lattice translations U x and 
U y must commute for all x,y 6 L. 

The product VU is a valid local, unitary quantum operation. The resulting global update rule is well defined and 
space homogeneous. The set of states includes a quiescent state as well as an "ancillary" set of states / subspace which 
can store the result of the "read" operation. The initial state of a local-unitary QCA consists of identical k d blocks of 
cells initialized in the same state. Local-unitary QCA are universal in the sense that for any arbitrary quantum circuit 
there is a local- unitary QCA which can simulate it. In addition any local- unitary QCA can be simulated efficiently 
using a family of quantum circuits (jl2n . Adding an additional memory register to each cell allows this class of QCA 
to model any reversible QCA of the Schumacher /Werner type discussed above. 



C. Block-partitioned and nonunitary QCA 

Brcnnen and Williams introduce a model of QCA which allows for unitary and nonunitary rules (fl7l ). 

Definition (Block-partitioned QCA). A Block-partitioned QCA is a 4-tuple {L,Y,,AT,M} consisting of (1) a 1- 
dimensional lattice of n cells indexed L = 0, .. . ,n — 1, (2) a 1-dimensional state space S, (3) a neighborhood scheme 
J\f, and (4) an update rule M applied over M . 
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Given a system with nearest-neighbor interactions, the simplest unitary QCA rule has radius r = 1 describing a 
unitary operator applied over a three-cell neighborhood j — j + 1: 

M (uoo,«oi, ttio, «n) = 100} <00| ® u o + |01) <01| ® u 01 + |10) (10| ® u 10 + |H) (H| ® «n , (4) 

where |a&) (o&| ® u a j, means update the qubit at site j with the unitary u a & if the qubit at the site j — 1 is in state 
\a) and the qubit at site j + 1 is in state \b). M commutes with its own 2-site translation. Thus, a partitioning 
is introduced by updating simultaneously all even qubits with rule M before updating all odd qubits with rule M. 
Periodic boundaries are assumed. However, by addressability of the end qubits simulation of a block-partitioned QCA 
by a QCA with boundaries can be achieved. 

Nonunitary update rules correspond to completely positive maps on the quantum states where the neighboring 
states act as the environment. Take a nearest-neighbor 1-dimensional Block-partitioned QCA. In the density operator 
formalism each quantum system p is given by the probability distribution p = J^iPi W) ("^l over outer products of 
quantum states \ip}- A completely positive map S(p) applied to state p is represented by a set of Krauss operators F lll 
which are positive operators that sum up to the identity F^F^ = 1. The map Sj b (p) acting on cell j conditioned 
on state a of the left neighbor and state b of the right neighbor can then be written as 

Sf G>) = |o6) {ab\ ® J2 FfpFf ® \ ab ) ( ab \ ■ ( 5 ) 

A* 

As an example, the CA rule '110' can now be translated into an update rule for cell j in a block-partitioned 
nonunitary QCA: 

F( = |00) (00| 8> V + |10) (10| <8> V + |11) (11| <g> a j x + |01) (01| ® |1)^ (1| (6) 
2^ = |01><01|® 1 1)^.<0| , (7) 
where a x is the Pauli operator. 

The implementation of such a block-partitioned nonunitary QCA is proposed in form of a lattice of even order 
constructed with an alternating array of two distinguishable species ABABABAB . . . that are globally addressable 
and interact via the Ising interaction. Update rules that generate and distribute entanglement are studied in this 
framework (fl7t). 

D. Continuous-time QCA 

Vollbrccht and Cirac initiate the study of continuous-time QCA (27). They show that the computability of the 
ground state energy of a translationally invariant n-neighbor Hamiltonian is QMA-hard. Their QCA model is taken 
up by Nagaj and Wocjam (f28h who use the term Hamiltonian QCA. 

Definition (Hamiltonian QCA). A Hamiltonian QCA is a tuple {(L, E = S p x Tig, )} consisting of (1) a 1- dimensional 
lattice of length L, (2) a finite set of orthogonal basis states S = H p x containing (2a) a data register E^, and (2b) 
a program register E p . 

The initial state encodes both the program and the data, stored in separate subspaces of the state space: 

L 

\4>) = (g) (\Pj) ® \dj))j (8) 
j'=i 

The computation is carried out autonomously. Nagaj and Wocjam show that, if the system is left alone for a 
period of time t = 0(L logL), polynomially in the length of the chain, the result of the computation is obtained 
with probability p > 5/6 — 0(l/log L). Hamiltonian QCA are computationally universal, more precisely they are in 
the complexity class BQP. Two constructions for Hamiltonian QCA are given in (|28l). one using a 10-dimensional 
state space, and the resulting system can be thought of as the diffusion of a system of free fermions. The second 
construction given uses a 20-dimensional state space and can be thought of as a quantum walk on a line. 

E. Examples of QCA 

Raussendorf presents an explicit construction of QCA and proves its computational universality (fl5l). The QCA 
lives on a torus with a 2 x 2 Margolus partitioning. The update rule is given by a single 4-qubit unitary acting on 
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2x2 blocks of qubits. The four-qubit unitary operation consists of swap operations, the Hadamard transformation, 
and a phase gate. The initial state of the QCA is prepared such that columns encode alternatingly data and program. 
When the QCA is running the data travel in one direction while the program (encoding classical information in 
orthogonal states) travels in the opposite direction. Where the two cross the computation is carried out through 
nearest-neighbor interaction. After a fixed number of steps the computation is done and the result can be read out of 
a dedicated "data" column. This QCA is computationally universal, more precisely, it is within a constant as efficient 
as a quantum logic network with local and nearest-neighbor gates. 

Shepherd, Franz, and Werner compare classically controlled QCA to autonomous QCA (fl6Th The former is con- 
trolled by a classical compiler that selects a sequence of operations acting on the QCA at each time step. The latter 
operates autonomously, performing the same unitary operation at each time step. The only step that is classically 
controlled is the measurement (and initialization). They show the computational equivalence of the two models. Their 
result implies that a particular quantum simulator may be as powerful as a general one. 

IX. COMPUTATIONALLY UNIVERSAL QCA 

Quite a few models have been shown to be computationally universal, i.e. they can simulate any quantum Turing 
machine and any quantum circuit efficiently. A Watrous-QCA simulates any quantum Turing machine with constant 
slowdown (0). The QCA defined by Van Dam is a finite version of a Watrous QCA and is computationally universal 
as well (@). Local- unitary QCA can simulate any quantum circuit and thus are computationally universal ([12). 
Block-partitioned QCA can simulate a quantum computer with linear overhead in time and space (|17). Continuous- 
time QCA are in complexity class BQP and thus computationally universal j27|). The explicit constructions of 
2-dimensional QCA by Raussendorf is computationally universal, more precisely, it is within a constant as efficient as 
a quantum logic network with local and nearest-neighbor gates (|15h . Shepherd, Franz, and Werner provide an explicit 
construction of a 12-state 1-dimensional QCA which is in complexity class BQP. It is universally programmable in 
the sense that it simulates any quantum-gate circuit with polynomial overhead (fl6l). Arrighi and Fargetton propose 
a 1-dimensional QCA capable of simulating any other 1-dimensional QCA with linear overhead (fl3l). 

Implementations of computationally universal QCA have been suggested by Lloyd (f29h and Benjamin (f30h. 

X. MODELING PHYSICAL SYSTEMS 

One of the goals in developing QCA is to create a useful modeling tool for physical systems. Physical systems that 
can be simulated with QCA include Ising and Heisenberg interaction spin chains, solid state NMR, and quantum lattice 
gases. Spin chains are perhaps the most obvious systems to model with QCA. The simple cases of such 1-dimensional 
lattices of spins are Hamiltonians which commute with their own lattice translations. Vollbrecht and Cirac have 
shown that the computability of the ground state energy of a translationally invariant rt-neighbor Hamiltonian is in 
complexity class QMA (H3). For simulating noncommuting Hamiltonians a block- wise update such as the Margolus 
partitioning has to be used (see Section lV.A|) . Here the fact is used that any Hamiltonian can be expressed as the sum 
of two Hamiltonians, H = H a + Hi,. H a and Hb can then, to a good approximation, be applied sequentially to yield 
the original Hamiltonian H, even if these do not commute. It has been shown that such 1-dimensional spin chains 
can be simulated efficiently on a classical computer (|3ll). It is not known, however, whether higher dimensional spin 
systems can be simulated efficiently classically. 

A. Quantum lattice gas automata 

Any numerical evolution of a discretized partial differential equation can be interpreted as the evolution of some CA, 
using the framework of lattice gas automata. In the continuous time and space limit such a CA mimics the behavior 
of the partial differential equation. In quantum mechanical lattice gas automata (QLGA) the continuous limit on a 
set of so called quantum lattice Boltzman equation recovers the Schrodinger equation j32Tj. The first formulation of a 
linear unitary CA is given in Ref. (f33h . Meyer coins the term quantum lattice gas automata (QLGA) and demonstrates 
the equivalence of a QLGA and the evolution of a set of quantum lattice Boltzman equations ((221 ; Meyer (|35h. 
Boghosian and Taylor |H), and Love and Boghosian |37t ) explored the idea of using QLGA as a model for simulating 
physical systems. Algorithms for implementing QLGA on a quantum computer have been presented in ((381: [39t l40l). 
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XI. IMPLEMENTATIONS 

A large effort is being made in many laboratories around the world to implement a model of a quantum computer. 
So far all of them are confined to a very finite number of elements and are no way near to a quantum Turing machine 
(which in itself is a purely theoretical construct but can be approximated by a very large number of computational 
elements). One existing experimental set-up that is very promising for quantum information processing and that does 
not suffer from this "finiteness" are optical lattices (for a review, see (41)). They possess a translation symmetry which 
makes QCA a very suitable framework in which to study their computational power. Optical lattices are artificial 
crystals of light and consist of hundreds of thousands of microtraps. One or more neutral atoms can be trapped in 
each of the potential minima. If the potential minima are deep enough any tunneling between the traps is suppressed 
and each site contains the same amount of atoms. A quantum register - here in form of a so-called Mott insulator 
- has been created. The biggest challenge at the moment is to find a way to address the registers individually to 
implement quantum gates. For a QCA all that is needed is implementing the unitary operation(s) acting on the 
entire lattice simultaneously. The internal structure of the QCA guarantees the locality of the operations. This is 
a huge simplification compared to individual manipulation of the registers. Optical lattices are created routinely 
by superimposing two or three orthogonal standing waves generated from laser beams of a certain frequency. They 
are used to study Fermionic and Bosonic quantum gases, nonlinear quantum dynamics, strongly correlated quantum 
phases, to name a few topics. 

A type of locally addressed architecture by global control was put forward by Lloyd (|29T). In this scheme a 1- 
dimensional array is built out of three atomic species, periodically arranged as ABC ABC ABC. Each species encodes a 
qubit and can be manipulated without affecting the other species. The operations on any species can be controlled by 
the states of the neighboring cells. The end-cells are used for readout, since they are the only individually addressable 
components. Lloyd showed that such a quantum architecture is universal. Benjamin investigated the minimum 
physical requirements for such a many-species implementation and found a similar architecture using only two types 
of species, again arranged periodically ABABAB (13fj|; I421 l43h . By giving explicit sequences of operations implementing 
one-qubit and two-qubit (CNOT) operations Benjamin shows computational universality. But the reduction in spin 
resources comes with an increase in logical encoding into four spin sites with a buffer space of at least four empty 
spin sites between each logical qubit. 

A continuation of this multi-species QCA architecture is found in the work by Twamley ([44]). Twamley constructs 
a proposal for a QCA architecture based on Fullerene (Ceo) molecules doped with atomic species 15 N and 31 P, 
respectively, arranged alternatingly in a one-dimensional array. Instead of electron spins which would be too sensitive 
to stray electric charges the quantum information is encoded in the nuclear spins. Twamley constructs sequences of 
pulses implementing Benjamin's scheme for one- and two-qubit operations. The weakest point of the proposal is the 
readout operation which is not well-defined. 

A different scheme for implementing QCA was suggested by Toth and Lent pH). Their scheme is based on the 
technique of quantum-dot CA. The term quantum-dot CA is usually used for CA implementations in quantum dots 
(for classical computation). The authors, therefore, call their model a coherent quantum-dot CA. They illustrate 
the usage of an array of N quantum dots as an N-qubit quantum register. However, the set-up and the allowed 
operations allow for individual control of each cell. This coherent quantum-dot CA is more a hybrid of a quantum 
circuit with individual qubit control and a QCA with constant nearest-neighbor interaction. The main property of a 
QCA, operating under global control only, is not taken advantage of. 

XII. FUTURE DIRECTIONS 

The field of QCA is developing rapidly. New definitions have appeared very recently. Since QCA are now considered 
to be one of the standard measurement-based models of quantum computation, further work on a consistent and 
sufficient definition of hig her-dimensional QCA is to be expected. One proposal for such a "final" definition has been 
put forward in Refs. (|14l ; \4w. 

In the search for a robust and easily implcmcntable quantum computational architectures QCA are of quite some 
interest. The main strength of QCA is global control without the need to address cells individually (with the possible 
exception of the read-out operation) . It has become clear that the global update of a QCA would be a way around 
practical issues related to the implementation of quantum registers and the difficulty of their individual manipulation. 

More concretely, QCA provide a natural framework for describing quantum dynamical evolution of optical lattices, 
a field in which the experimental physics community has made huge progress in the last decade. 

The main focus so far has been on reversible QCA. Irreversible QCA are closely related to measurement-based 
computation and remain to be explored further. 
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XIII. SECONDARY REFERENCES 

Summaries of the topic of QCA can be found in Chapter 4.3 of Gruska (j47t ). and in Rcfs. (|48l ; l49h. 
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